Back

The Journal of Physical Chemistry B

American Chemical Society (ACS)

Preprints posted in the last 7 days, ranked by how well they match The Journal of Physical Chemistry B's content profile, based on 167 papers previously published here. The average preprint has a 0.09% match score for this journal, so anything above that is already an above-average fit.

1
A Simple Method to Distinguish Active and Inactive Aptamers by Analyzing the Ruggedness of the Aptamer Free Energy Landscape

Subramanian, G.; Thiel, W.; Singh, R.

2026-08-29 bioinformatics 10.64898/2026.08.26.747184 medRxiv
Top 0.7%
3.2%
Show abstract

Aptamers are structured nucleic acid ligands capable of high affinity, high specificity molecular recognition generated using variations of the SELEX (Systematic Evolution of Ligands by Exponential Enrichment) process. However, SELEX often produces sequences that enrich yet may lack binding efficacy. We propose a measure called the Ruggedness Composite Index (RCI) along with a method for computing it, that can be used to distinguish binding-competent ('active') aptamers from weak or non-binding ('inactive') aptamers. Given a set of aptamers, RCI incorporates information on their fragmentation (landscape partitioning), basin entropy (metastable state distribution), cumulative density irregularity (non-uniform occupancy), and structural energy correlation length (structure-energy coupling scale). We test whether secondary-structure folding energy landscape topology distinguishes active from inactive aptamers using a multiscale level set framework across six datasets. Active aptamers show lower RCI values and occupy smoother, funnel-like conformational spaces, while inactive aptamers show higher RCI values, reflecting fragmented, high-entropy landscapes. By contrast, classical thermodynamic features, such as minimum free energy, show limited discrimination between active and inactive aptamers. In all datasets, sequences that exhibit enrichment which is not monotonic but lack specificity exhibit elevated ruggedness, indicating landscape topology can predict non-specific enrichment. These results indicate that folding landscape organization can be used as a predictor of aptamer activity and establish RCI as a simple, mechanistically interpretable measure for improving candidate prioritization, especially in therapeutic aptamer discovery.

2
Effects of Cholesterol on Nanodisc Formation and Magnetic Alignment in DMPC and Glycyrrhizic Acid Systems Probed by 31P and 14N Solid-State NMR

Rokonujjaman, M.; Wi, S.; Ramamoorthy, A.

2026-08-29 biophysics 10.64898/2026.08.26.747314 medRxiv
Top 0.8%
2.8%
Show abstract

Nanodiscs and bicelles are widely used as membrane mimetics for structural studies of membrane-associated systems. Studies have reported that their magnetic alignment behavior and phase stability are highly sensitive to composition and temperature. In this study, we systematically investigate the effects of cholesterol on bicelle formation and magnetic alignment in DMPC + 0.2 glycyrrhizic acid (GA) systems using a combined 31P and 14N solid-state NMR experimental and simulation-based approach. Temperature dependent 31P NMR spectra reveal a clear transition from vesicle dominant to aligned bicelles/nanodsics phase, while 1N quadrupolar splitting and lineshape analysis provides quantitative insights into heterogeneous lipid bilayer populations, distinguishing large aligned nanodiscs (B(L)), small nanodiscs (B(S)), and isotropic/random components (B(R)). A strong correlation is observed between the 31P derived bicelle fraction and the 14N B(L) population, confirming that macroscopic alignment in the presence of an external magnetic field directly reflects the growth of large, well-ordered nanodiscs. Cholesterol is found to play a critical dual role by modulating membrane order and curvature. At low cholesterol concentration (0 to 5 mole percent), nanodiscs alignment occurs gradually with increasing temperature, while at higher cholesterol concentration (15 to 25 mole percent), the alignment is delayed and accompanied by broader spectral features, indicating structural heterogeneity. Notably, 10 mole percent cholesterol consistently provides the optimal balance, enabling efficient temperature dependent conversion to aligned bicelles while maintaining high B(L) populations (about 70-80 percentage) and minimal isotropic fractions. In contrast, higher cholesterol maintains significant B(S) and B(R) populations, even at elevated temperature. The 14N quadrupolar coupling (Cq is approximately 8.5 to 9.2 kHz for aligned nanodiscs) remains nearly invariant across compositions, showing that cholesterol does not change local headgroup dynamics but instead redistributes lipid populations. These findings establish a combined 31P and 14N solid -state NMR approach provides a valuable platform for quantitatively correlating membrane structure, dynamics, and alignment, offering practical guidelines for optimizing bicelle systems for high resolution solid-state NMR studies of membrane associated biomolecules.

3
Towards transferable explicit-solvent coarse-grained models for biomolecular condensates

Toplek, F. B.; Borges-Araujo, L.; Lindorff-Larsen, K.; Everaers, R.; Souza, P. C. T.; Morozova, T. I.

2026-08-29 biophysics 10.64898/2026.08.27.747511 medRxiv
Top 0.8%
2.7%
Show abstract

Biomolecular condensates formed by intrinsically disordered proteins require molecular models that accurately describe proteins in both dilute solution and condensed phases. Explicit-solvent coarse-grained models offer an attractive balance between chemical resolution and computational efficiency. Yet, it remains unclear whether improving dilute-state properties is sufficient to obtain an accurate description of condensates. Here, we address this question by introducing minimal modifications to the Martini 3 force field that combine recent advances in bonded interactions with refined protein-water interactions and strengthened glycine self-interactions, while preserving the underlying chemical transferability of the model. The resulting model substantially improves the description of single-chain conformations across a diverse benchmark of disordered proteins. We then investigate phase separation of the well-characterized low-complexity domain of heterogeneous nuclear ribonucleoprotein A1 and its sequence variants. The model reproduces several key physicochemical properties of biomolecular condensates, including chain expansion in the dense phase, sequence-dependent intermolecular contacts, protein diffusion and its relation to single-chain dimensions, and hydration, while revealing quantitative limitations in condensate density, phase equilibria, and ion partitioning. Our results show that improving dilute-state behaviour translates into a better description of condensed-phase properties, including condensate density, but is not sufficient to quantitatively reproduce the equilibrium between the dilute and dense phases.

4
A thermodynamic framework for mapping elastic recoil mechanism across the human proteome

Desai, R.; Pople, D.; Musale, A.; Jain, S.; Sajjad, I.; Wittebort, R. J.; Koder, R. L.; Nanda, V.

2026-08-30 biophysics 10.64898/2026.08.28.747957 medRxiv
Top 1%
1.9%
Show abstract

The folding thermodynamics of proteins are dominated by two opposing forces, the loss in backbone entropy and the packing of hydrophobic groups. The same forces are major contributors to the extension thermodynamics of elastic proteins with the distinction that both processes act in concert, favoring the higher chain and solvent entropy of a relaxed conformation. The relative entropic contributions specify the recoil mechanism; human elastin recoil is primarily driven by hydrophobic forces, whereas fly resilin has a rubber-like mechanism driven by backbone entropy. Despite the importance of elastic proteins to tissue biomechanics, few have been identified, let alone characterized to the same extent as elastin and resilin. We develop a thermodynamic framework that maps proteins by sequence-derived estimates of extension-induced backbone and solvent entropy changes. Putative elastic proteins are proposed and classified by recoil mechanism based on estimated thermodynamic features. Proteins that map to elastic regions are overrepresented by the skin proteome. The set of predicted elastic domains is further extended by incorporating sequence context embedded in protein language models. Protein domains with distinct thermodynamic recoil mechanisms cluster on the latent space manifold. Some of these domains are anticipated to have roles within molecular machines, expanding the scope of elastic protein function beyond mechanical materials like elastin and resilin.

5
Geometric characterization of the HSV - 1 glycoprotein B - amyloid β interaction in Alzheimer's disease using Forman-Ricci curvature

Bou Dagher, L.; Han, Z.; Zhou, S.; Fülöp, T.; Desroches, M.; Rodrigues, S.

2026-08-29 bioinformatics 10.64898/2026.08.26.747308 medRxiv
Top 1%
1.5%
Show abstract

Alzheimer's disease is characterized by the accumulation and aggregation of amyloid-{beta}(A{beta}), but the molecular mechanisms linking environmental and infectious factors to A$\beta$ conformational changes remain incompletely understood. Herpes simplex virus type 1 (HSV-1) has been proposed as a potential contributor to AD pathology, and interactions between the viral glycoprotein B (gB) and A$\beta$ may influence the conformational behaviour of the peptide. Molecular dynamics (MD) simulations provide atomic-scale information on such interactions, but conventional structural descriptors may not fully capture changes in the organization of residue interaction networks. Here, we introduce a graph-geometric framework based on Forman-Ricci curvature to characterize the evolution of residue interaction networks during MD simulations. Each simulation frame is represented as a residue interaction graph based on C--C contacts, and residue-wise curvature profiles are analysed across time. We apply the framework to A{beta}1-42 in isolation and in complex with HSV-1 gB. Conventional MD analyses indicate stable association of the simulated complex, favourable interaction energetics, and conformational changes in A{beta}, including a transition from -helical structure toward {beta}-turn-rich conformations over the simulated timescale. Forman-Ricci curvature reveals pronounced and spatially localized remodelling of the A{beta} residue interaction network in the complex, with the strongest changes concentrated in the C-terminal region. These regions also exhibit reduced temporal curvature fluctuations and progressively distinct geometric behaviour throughout the simulation. Hierarchical clustering further identifies cooperative groups of residues with coordinated curvature dynamics, including a prominent C-terminal domain. Together, these results demonstrate that Forman-Ricci curvature provides a complementary description of biomolecular dynamics by capturing changes in the geometric organization of residue interaction networks that are not directly represented by conventional structural descriptors. The framework provides a general computational approach for studying network-level structural remodelling in protein molecular dynamics and offers a quantitative perspective on the conformational consequences of HSV-1 gB--A{beta} association.

6
Heterotypic interactions and sequence features modulate cellular reflectin condensate dynamics

Phan, C.; Watanabe, R.; Le, V. Q.; Walsh, S.; Levenson, R.

2026-08-29 biophysics 10.64898/2026.08.27.747642 medRxiv
Top 2%
1.0%
Show abstract

Reflectin proteins drive dynamic structural coloration in cephalopods by organizing into dense intracellular lamellar structures that dictate local refractive index. While reconstituted reflectins readily undergo liquid-liquid phase separation in vitro, these assemblies frequently undergo dynamic arrest, vitrifying into non-dynamic condensates. Here, we investigate the primary sequence features, post-translational modifications, and heterotypic interactions that regulate the material properties of reflectin condensates within the crowded cellular environment of mammalian HeLa cells. Using confocal microscopy and fluorescence recovery after photobleaching (FRAP), we demonstrate that canonical block copolymeric A-type reflectins readily form dynamically arrested condensates, with the linker blocks primarily responsible for the observed arrest. In contrast, non-canonical B/C reflectin variants exhibit significantly greater fluidity and rapid recovery kinetics. We show that phosphomimetic substitutions progressively fluidize some reflectin condensates. Lastly, we find that heterotypic condensates composed of canonical and non-canonical reflectins in combinations associated with reversible iridescence in squid substantially enhance canonical mobility. Our findings establish a biophysical framework in which phosphorylation and heterotypic mixing cooperatively suppress dynamic arrest, enabling the reversible material transitions required for active cephalopod camouflage and communication.

7
Amyloid Polymorphism of Lysozyme Governs Cross-Seeding of Insulin Aggregation

Metkar, S.; Eerati, V.; Ramamoorthy, A.

2026-08-30 biophysics 10.64898/2026.08.26.747312 medRxiv
Top 2%
0.5%
Show abstract

Amyloid fibrils are highly ordered protein aggregates characterized by a conserved cross-{beta}-sheet architecture despite originating from structurally diverse precursor proteins. Growing evidence suggests that interactions between different amyloidogenic proteins can modulate aggregation pathways through heterologous cross-seeding; however, the influence of seed polymorphism on the structure and biological properties of cross-seeded fibrils remains poorly understood. Here, we investigated the cross-seeding of native human insulin by two structurally distinct polymorphs of hen egg-white lysozyme (HEWL): flexible fibrils (FFs) and rigid fibrils (RFs). Native insulin remained stable under physiological conditions and underwent spontaneous fibrillation only under acidic conditions. In contrast, both HEWL polymorphs efficiently induced insulin aggregation at physiological pH, bypassing the nucleation barrier. Thioflavin T fluorescence, circular dichroism spectroscopy, and transmission electron microscopy revealed that lysozyme FFs templated the formation of insulin flexible fibrils (IFFs), whereas lysozyme RFs produced insulin rigid fibrils (IRFs), demonstrating that the structural characteristics of the parental HEWL polymorphs were propagated during heterologous cross-seeding. The toxicity of the resulting insulin fibrils was evaluated in SH-SY5Y neuronal cells and CCF-STTG1 astrocytes. IFFs exhibited minimal cytotoxicity and only subtle morphological alterations, whereas IRFs caused modest reductions in cell viability accompanied by more pronounced cellular damage. These findings demonstrate that the structural polymorphism of HEWL fibrils governs both the architecture and biological activity of cross-seeded insulin fibrils, highlighting amyloid polymorphism as an important determinant of heterologous amyloid propagation and a potential design principle for engineering functional amyloid-based biomaterials and protein delivery platforms.

8
Engineering a highly active thermophilic F1-ATPase by homolog-guided exploration and machine-learning-assisted prioritization

Kobayashi, R.; Miyake, K.; Oya, T.; Ueno, H.; Saito, Y.; Noji, H.

2026-08-29 biophysics 10.64898/2026.08.27.747693 medRxiv
Top 2%
0.5%
Show abstract

The rotary motor F1-ATPase has been extensively studied as a model molecular machine, yet rational engineering of its catalytic activity remains challenging because ATP hydrolysis is regulated by long-range intersubunit allostery and large conformational transitions. Here, we developed a homolog-guided engineering strategy to increase the maximum rotation rate of the thermophilic Bacillus PS3 F1-ATPase (TF1). Candidate mutation sites were first identified by comparing TF1 with the homologous enzymes bovine mitochondrial F1 (bMF1) and Paracoccus denitrificans F1 (PdF1), both of which exhibit higher maximum rotation rates than TF1. Systematic exploration of these sites identified four activity-enhancing hotspots, followed by focused hotspot exploration and machine-learning-assisted prioritization of combinatorial mutants. The best mutant, TF1({beta}Y313L/{beta}E332S), exhibited a 1.8-fold higher maximum rotation rate than TF1(WT) while retaining its functional thermostability. Interestingly, activity-enhancing substitutions were not limited to the residues conserved in both bMF1 and PdF1, indicating that the bMF1-PdF1 consensus substitutions effectively identify activity-enhancing hotspots rather than uniquely defining the optimal amino acid. Machine-learning-assisted exploration efficiently prioritized highly active mutants, although the predictive performance was limited by the relatively small training dataset and epistatic interactions among mutations. Kinetic and structural comparisons further provided mechanistic insights into the enhanced catalytic activity of the engineered mutant. Together, these results establish a practical strategy for engineering complex molecular motors by combining homolog-guided hotspot identification with focused hotspot exploration.

9
Modeling Dynamics of Contact Inhibition of Proliferation and Structural Order in a Confluent Epithelium

Ghosh, J.; Bhattacharjee, T.; Dutta, S.

2026-08-29 biophysics 10.64898/2026.08.26.747344 medRxiv
Top 2%
0.5%
Show abstract

Contact inhibition of proliferation (CIP) enables epithelial tissues to self-regulate growth and maintain tissue homeostasis. However, how cell-level mechanical contact, tissue-scale structural order, and proliferation kinetics interplay remains a fundamental open question in living matter physics. Here, we present a particle-based model of a confluent epithelial monolayer governed by overdamped dynamics, where individual cells interact via a two-dimensional hard core- soft shoulder potential. By comparing structural evolution during quasistatic densification with previously reported experimental division kinetics, we find that the dynamics of proliferation arrest mimics the onset of direct steric contacts between the hard cores of the shell. Identifying hard core contacts as the physical driver of CIP, we couple our mechanical model with a stochastic Monte Carlo division scheme in which the instantaneous division rate decreases to zero from an intrinsic value as the number of hard core contact increases to six from zero. We demonstrate that for high intrinsic division rates, the cellular densification outpaces mechanical relaxation. This kinetic mismatch drives premature hard-core contact formation, shifts the onset of jamming and contact inhibition to lower packing fractions, and induces increasingly disordered transient configurations before the tissue universally converges to a hexagonal close-packed limit. Our model's predicted division kinetics and structural order evolution are consistent with epithelial monolayer experiments, both reported and our own. This minimal physical framework links single-cell steric contact mechanics directly to tissue-scale growth regulation and structural evolution.

10
Single-Molecule Nanopore Profiling of p53-TAD Conformational Dynamics, Interactions, and Inhibition

DeCoeur, D.; Schultz, S.; Chen, J.; Chen, M.

2026-08-29 biophysics 10.64898/2026.08.28.747917 medRxiv
Top 3%
0.3%
Show abstract

Investigating the conformational dynamics of intrinsically disordered proteins (IDPs) is essential to understanding how their structural heterogeneity underlies function and how their dysregulation contributes to diseases. Here, we utilized an MspA nanopore-based approach for studying the conformational dynamics and interactions of IDPs at the single-molecule level. The platform was demonstrated using the intrinsically disordered transactivation domain of tumor suppressor p53 (p53-TAD), one of the important proteins in cancer biology. We showed that MspA can stably capture p53-TAD and resolve up to six distinct current states with frequent interconversions, revealing a rich conformational landscape. The nanopore also detected the effect of a cancer-associated double mutational variant, N29K/N30D. Combining experiments with steered molecular dynamics simulations, we showed that the mutant sampled compact conformational states more frequently than wild type, consistent with previous NMR studies. Importantly, the MspA platform enabled direct monitoring of E3 ligase MDM2 binding to p53-TAD and resolved how this interaction is inhibited by anti-cancer compound epigallocatechin gallate (EGCG). Notably, EGCG stabilizes one of the six states sampled by p53-TAD, providing a mechanistic explanation for its inhibitory effect. Together, these findings demonstrate the promise of the nanopore platform for label-free monitoring of IDP conformational dynamics, modulation, binding and inhibition at single-molecule resolution.

11
Prot2Surf: fast analysis of protein - surface binding modes

Muniz-Chicharro, A.; Tanriver, G.; Gora, A.

2026-08-29 bioinformatics 10.64898/2026.08.26.747352 medRxiv
Top 3%
0.3%
Show abstract

Summary: Prot2Surf is a software tool designed for the characterization and prediction of protein association to surfaces. In this application note, Prot2Surf was tested using catalytic domains of the lytic polysaccharide monooxygenases (LPMOs), interacting with native surfaces. The results show that the software can efficiently analyze key binding features, including protein-surface distances, distances between catalytically reactive atoms, and the orientation angle between surface chains and the protein. These features are essential for distinguishing productive binding poses in these protein-surface systems and for understanding interaction patterns that provide guidance on protein engineering. Prot2Surf performs these analyses within seconds to a few minutes, providing a fast and accessible framework to post-process and characterize protein-surface encounter complexes. Availability and implementation: Prot2Surf, which is written in Fortran90, is documented and freely available as open source on GitHub: https://github.com/TUNNELING-GROUP/Prot2Surf. In order to run Prot2Surf, users should also install the SDA software package which is freely available at https://www.h-its.org/downloads/sda7/.

12
A time-delayed mechanochemical feedback model reconciles stable maintenance and dynamic remodeling of cell-matrix adhesions

Matsumoto, E.; Yokoyama, S.; Matsui, T. S.; Araki, T.; Deguchi, S.

2026-08-30 biophysics 10.64898/2026.08.28.747716 medRxiv
Top 3%
0.3%
Show abstract

Focal adhesions maintain force-bearing attachment between cells and the extracellular matrix but can also undergo dynamic remodeling. Their assembly and actomyosin tension are coupled through mechanochemical feedback. The processes underlying this feedback are not instantaneous and therefore involve a time delay. However, how this delayed feedback gives rise to stable adhesion maintenance or dynamic remodeling remains unclear. Here, paired time-lapse measurements of vinculin fluorescence and traction stress revealed distinct local adhesion-force dynamics, including low-fluctuation and recurrent fluctuation patterns. To examine how these patterns could arise, we formulated a minimal mechanochemical model coupling focal adhesion assembly and actomyosin force through delayed reciprocal feedback. The model exhibited stable and oscillatory modes depending on feedback strength, the balance of opposing feedback effects, and the effective feedback delay. Bistability and hysteretic switching also occurred in a subset of parameter space, and the oscillation period followed a power-law relation with the delay. These results suggest that stable adhesion maintenance and dynamic remodeling can emerge from a common mechanochemical feedback architecture.

13
FlexiTAC enables controllable PROTAC linker generation across diverse structural settings using a Bayesian flow network with posterior guidance

Li, Y.; Zhao, Y.; Zhou, L.; Huang, C.; Xu, Q.; Chen, Y.; Qin, Z.; Fan, K.; Yang, J.; Cao, D.

2026-08-30 bioinformatics 10.64898/2026.08.26.747172 medRxiv
Top 3%
0.3%
Show abstract

Linker chemistry and conformation are central determinants of PROTAC activity, shaping ternary-complex geometry, cooperativity, target-lysine presentation and cellular permeability. Existing linker generators often lack explicit control over linker flexibility, require predefined attachment sites and linker lengths, or produce structures that demand substantial geometric correction, limiting their utility in practical PROTAC design. Here we introduce FlexiTAC, a Bayesian flow network that jointly generates linker atom types and coordinates from the warhead and E3-ligase-ligand contexts. We also assemble PROTAC-3D, a quality-controlled collection of 63,554 component-resolved PROTAC structures for model training, and PROTAC-Bench, which covers molecular quality, fragment preservation, geometric fidelity, conformational stability, fragment awareness, rediscovery and sampling efficiency. Compared to the best 3D baseline models, FlexiTAC improves validity by 12.0-12.7% and achieves the highest PoseBusters pass rate of 79.5%-80.0%. A differentiable guidance module shifted generated linkers along a conformational ensemble-derived rigidity axis without retraining the generator. In silico case studies further show that the model can accept crystal-derived, redocked or predicted structural inputs. Together, FlexiTAC, PROTAC-3D and PROTAC-Bench establish an integrated and reproducible framework for data-driven PROTAC linker design, combining controllable structure-conditioned generation with standardized training data and evaluation protocols. This framework expands the linker chemical and conformational space accessible to computational exploration, provides a foundation for future method development and enables the systematic generation of structure-conditioned linker designs with tunable conformational flexibility.

14
Biochemical and Binding Characterization of a Riboflavin Analogue Tethered to Biotin

Marincean, S.; Smith, S. R.; Branscum, T.; Ratajczak, A.; Benore, M. A.

2026-08-31 biochemistry 10.64898/2026.08.29.748002 medRxiv
Top 4%
0.1%
Show abstract

The binding affinities of a chimeric analog of a riboflavin derivative linked to biotin, (6- (7,8-dimethyl-2,4-dioxo-3,4-dihydrobenzo[g]pteridin-10(2H)-yl)hexyl 5-((3aS,4S,6aR)-2- oxohexahydro-1H-thieno[3,4-d]imidazol-4-yl)pentanoate), referred to as C6-Rf-biotin-tag, to the riboflavin binding retain or streptavidin are in the M range, 1.29 {+/-} 0.277 and 3.00 {+/-} 0.459, respectively. These values suggest that C6-Rf-biotin-tag has potential applications in diagnostic assay and labelling target flavin binding proteins. The C6-Rf-biotin-tag which was characterized with respect to physical and biochemical properties retains UV/Vis spectroscopic and fluorescence behavior similar to riboflavin.

15
Cryo-EM structure of CYP2C9 reveals a dimer-of-trimers assembly

Tanino, H.; Tsujino, H.; Nakao, T.; Oie, C.; Makino, F.; Miyata, T.; Kasai, K.; Namba, K.; Inoue, T.

2026-08-31 biophysics 10.64898/2026.08.30.747968 medRxiv
Top 4%
0.1%
Show abstract

Human cytochrome P450 2C9 (CYP2C9) is a hepatic microsomal enzyme involved in the oxidative metabolism of clinically important drugs, but the structural organization of its oligomeric assemblies outside crystallographic packing environments remains poorly understood. Here, we report the cryo-EM structure of human CYP2C9 determined under aqueous, membrane-free conditions at 3.31 Angstrom resolution. The structure reveals a C2-symmetric hexameric assembly organized as a dimer of trimers. Individual protomers retain the conserved P450 fold and heme-binding architecture observed in previously reported crystal structures, indicating that assembly formation does not substantially perturb the catalytic core. The hexamer is stabilized by defined intra-trimer interfaces involving the N-terminal region and residues around Trp212 and Phe482, together with inter-trimer interfaces involving Leu71 and the 220-227 loop. These interfaces are distinct from the crystal packing contacts observed in CYP2C9 crystal structures, demonstrating that the assembly is not a simple recapitulation of crystallographic packing. Notably, the inter-trimer interface is located near the FG-loop-containing surface previously implicated in membrane association. This suggests that the observed hexamer may represent a membrane-free association of two trimers through membrane-related surfaces, whereas the trimeric arrangement itself may be compatible with membrane-associated organization. The structure therefore provides a framework for investigating how trimer formation, membrane interaction and local conformational changes in the FG-loop region may influence CYP2C9 function.

16
Mechanism of heme binding by CP motifs in the BACH1 DNA-binding region

Huang, Y.; Fairall, L.; Muskett, F. W.; Dominguez, C.; Hudson, A.; Schwabe, J. W.

2026-08-31 biochemistry 10.64898/2026.08.28.747782 medRxiv
Top 4%
0.1%
Show abstract

BACH1 is a heme-regulated basic-leucine-zipper containing transcriptional repressor that binds its DNA recognition elements as a heterodimer with MAFK. Heme-binding is thought to be mediated by several Cys-Proline (CP) motifs and this results in dissociation of the heterodimer from DNA. The mechanism of heme-binding and heme-mediated DNA dissociation remains unresolved. We have used UV-visible spectroscopy, 2D-NMR and DNA-binding assays to explore both heme-binding and DNA dissociation of a minimal BACH1 construct containing 2 CP motifs (C492(CP5) and C646(CP6)) flanking the DNA-binding domain. We find that heme is able to bind to both CP motifs, but also to other non-CP cysteines and histidines in the construct. Using NMR spectroscopy, we identify a structured binding pocket in which heme interacts with both C646(CP6) and Cys621. However, DNA-binding assays show that C646(CP6) is not required for heme-mediated DNA dissociation of the BACH1:MAFK heterodimer. Using UV-visible spectroscopy we show that C492(CP5) also recruits heme with a second ligand, a conserved histidine, His559, in the BACH1 DNA-recognition helix. Mutation of C492(CP5) reduces but does not abolish heme-mediated dissociation from DNA. Our findings suggest a mechanism for heme-binding to BACH1 and heme-mediated dissociation from DNA.

17
OmniScore: Universal Scoring of Diverse Biomolecular Complexes via Equivariant Geometry-Aware Discrete Representation Learning

Bui, T.-C.; Lee, J.; Ko, J.

2026-08-29 bioinformatics 10.64898/2026.08.28.747942 medRxiv
Top 4%
0.1%
Show abstract

Scoring biomolecular complexes is central to structure assessment and drug discovery, yet the complexes themselves vary widely in pose, size, and molecular composition. A scoring function tuned for one interaction type rarely carries over to another, and most existing methods compound the problem by leaning heavily on task-specific labels. We introduce OmniScore, a universal structure-based framework that learns a shared geometry-aware representation of complexes once and then adapts it to downstream scoring through lightweight task-specific heads. OmniScore couples a graph view and a sequence view of each structure, encodes its three-dimensional geometry, and compresses representations into a compact latent space that a reconstruction module and prediction heads can reuse. We pretrain this backbone on diverse datasets including complexes, monomers, and small molecules with complementary objectives: coordinate recovery, correcting corrupted input tokens, predicting molecular identity, and grounding the representation in structure-level physical quantities. Across the evaluated benchmarks, OmniScore gave the best antibody-antigen and nanobody-antigen quality assessment on all reported metrics compared to state-of-the-art baselines. Its frozen residue embeddings matched the state-of-the-art protein-tokenization method with an average functional-site accuracy of 71.8% on a standard residue-level benchmark. On protein-ligand scoring and ranking benchmarks, it performed on par with methods built specifically for that single task. These results suggest that geometry-aware pretraining can provide a reusable scoring backbone for tasks that depend on interfacial and residue-level structure, within the evaluated settings.

18
Chemi-Proteome Language Attention Network Empowers Fragment-Based Ligand Interactome and Binding Sites Discovery with Evidence

Liao, B.; He, J.; zhao, M.; Cui, X.; Cui, Y.; Dong, C.; Sun, H.; Zhang, L.; Zhang, J.

2026-08-30 bioinformatics 10.64898/2026.08.26.747036 medRxiv
Top 4%
0.1%
Show abstract

Deep learning has accelerated drug discovery, yet most existing models are trained using in vitro affinity datasets and consequently remain disconnected from the cellular context in which functional ligand-protein interactions occur. This limitation hinders the ability to reflect the complexity of native interactomes and characterize biological responses to molecular perturbation. Here we introduce C-PLANK (Chemi-Proteome Language Attention NetworK), a deep learning framework trained on fragment-protein interactions profiled directly in living cells using fully functionalized fragment (FFF) chemoproteomics. C-PLANK combines physicochemical embeddings with a bilinear attention network (BAN) to model both global cellular context and local residue-atom interactions, generating interpretable interaction fingerprints. Particularly, C-PLANK incorporates Cellular Interaction State Index (CISI), a systems-level evidential metric that contextualizes the biological plausibility of each predicted interaction against the global cellular interaction landscape. Across 431 ligand interactomes curated from eight independent chemoproteomic studies, C-PLANK consistently outperformed current state-of-the-art interaction prediction frameworks under both random and cold-protein evaluation settings. The inferred interaction fingerprints aligned with orthogonal evidence from structure-based pocket predictions, co-crystal structures, and cellular binding-site annotations. C-PLANK further generalized to unseen ligands. In a cellular target-focused discovery campaign, C-PLANK identified a previously unrecognized ligand that was subsequently advanced into an active chemical probe acting as a SIRT3 agonist in cellular assays. By learning directly from cellular chemoproteomics, C-PLANK moves beyond isolated interaction prediction toward cellular interaction-state modelling, establishing a computational foundation for future digital-twin frameworks in drug discovery.

19
Cereblon on Steroids: Beyond the Canonical Ligand Space

Herrmann, A.; Heim, C.; Maiwald, S.; Boichenko, I.; Neuenschwander, M.; Oder, A.; Hernandez Alvarez, B.; Lupas, A. N.; von Kries, J. P.; Hartmann, M. D.

2026-08-31 biochemistry 10.64898/2026.08.28.747849 medRxiv
Top 5%
0.1%
Show abstract

Cereblon (CRBN) is widely used in targeted protein degradation, but its ligand space has remained dominated by a narrow set of cyclic imide chemotypes. Here, we show that the accessible CRBN ligand space extends substantially beyond this canonical space. A high-throughput screen of > 40,000 compounds, followed by orthogonal biophysical validation, X-ray crystallography and SAR analyses, identified several chemically distinct ligand classes. These include linear acetyl-based motifs, a phthalide-derived scaffold, steroidal compounds, and a range of bicyclic ligands. They engage CRBN through distinct recognition modes, several of which deviate from the canonical hydrogen-bonding pattern. Steroidal scaffolds were particularly notable: cortisone binds the human CRBN thalidomide-binding domain with an affinity comparable to thalidomide, with its A-ring occupying the tri-tryptophan pocket in a glutarimide-like orientation despite lacking the canonical imide NH donor. SAR within this series showed substantial tolerance for chemical modification and scaffold simplification, raising the possibility that endogenous steroidal metabolites may contribute to the physiological ligand landscape of CRBN. Bicyclic lactams additionally provided synthetically accessible scaffolds with tunable affinity and promising sites for linker attachment. Across the identified ligand classes, none of the tested representatives induced detectable degradation of canonical CRBN neosubstrates, and several showed largely clean proteomic profiles. Together, these findings broaden the chemical, mechanistic and potential physiological landscape of CRBN recognition and provide diverse starting points for alternative, potentially neosubstrate-sparing CRBN recruiters.

20
Subtyping active-site inhibitor binding mode to Abl kinase using super-resolution nanopore tweezers.

Ly, N.; Wang, Y.-H.; Foster, J.; DeCoeur, D.; Nguyen, L.; Wu, B.; Milenkovic, O.; Chen, M.

2026-08-29 biophysics 10.64898/2026.08.27.747610 medRxiv
Top 5%
0.1%
Show abstract

Accurate determination of kinase inhibitor binding modes could provide essential information for understanding resistance mechanisms and accelerating drug discovery. While conventional structural methods such as X-ray crystallography, cryo-EM and NMR provide high-resolution information but are low-throughput and capture largely static snapshots of dynamic protein-ligand interactions Here, we introduce a single-molecule nanopore tweezer platform that functionally subtypes ATP-competitive Abl kinase inhibitors by resolving distinct ionic current signatures of Abl-inhibitor complexes. This approach distinguishes Type I, Type IIA, and Type IIB inhibitors without structural determination. We further show how clinically relevant Abl variants (T315I and E255V) reshape inhibitor engagement and binding modes. By combining baseline probability features with wavelet-based time-frequency descriptors, ensemble machine-learning models achieved 97.5% classification accuracy across seven kinase inhibitor binding modes at sub-angstrom resolution and enabled deconvolution of mixed-inhibitor samples at nanomolar concentrations. These results establish nanopore tweezers as a label-free, super-resolution platform for profiling kinase conformational states and inhibitor binding modes, complementing structural approaches and supporting precision oncology.